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Abstract 

The paper establishes tight lower bound for effective elastic stress energy of 
two-dimensional three-phase anisotropic composites. It is assumed that the 
materials are mixed with fixed volume fractions and that one of the phases is 
degenerated to void, i.e. the effective material is porous. The bound expands 
the Hashin-Shtrikman and translation bounds to multiphase structures. It 
is derived using a combination of the translation method and additional in- 
equalities on the stress fields in materials; similar technique was used by 
Cherkaev (Mech. Mater. 41, 411-433, 2009) for isotropic multiphase com- 
posites. It is shown that in the high-porosity regime the obtained energy 
estimate is exact and it is realized on microstructures in the form of high- 
rank laminates whose optimal geometrical parameters are explicitly found. 
To this end, the field matching method is elaborated for the plane elasticity 
problem, following Albin et al. (J. Mech. Phys. Sohds 55, 1513-1553, 2007) 
and Cherkaev and Zhang (Int. J. Solids Struct. 48, 2800-2813, 2011) where 
similar technique was developed for conductivity case. Corresponding for- 
mulae solving the "G-closure problem" of a three-phase composite are also 
derived. Conjectures regarding low-porosity regimes are presented. However, 
full discussion of this issue is postponed to a separate publication. 
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1. Introduction 



Significance of the problem Optimization of composite microstruc- 
tures is important today because teclmological capabilities allow for manufac- 
turing a huge variety of microscopic designs for roughly the same price, and 
one wants to know what "the best" microstructure is. There is no boundary 
between optimal structural design in classical engineering sense and optimal 
composite material as the latter is also a structure at microlevel. Optimal 
large-scale structures are made from optimal microstructures (composites) 
and the main difference between them is that the composite problem is solved 
for a periodic domain and periodic boundary conditions, which permits for 
an explicit solution. Besides the optimal structures, one wants to know the 
range of the improvement of overall composite properties by varyi ng the mi- 
crostructure - the quas iconvex envelope problem, see for example (jCherkaevl . 



2000l : iDacorognaL 120081 ) - which opens ways to construction of metamaterials. 



i.e. structures with unusual responses. 

So far, the vast majority of available results deals with two-material com- 
posites. Meanwhile, numerous applications call for optimal design of multi- 
material composites, or even porous composites from two elastic materials 
and void, especially applications that utilize multi-physics, i.e. elastic and 
electromagnetic properties and those that deal with structures best adapted 
to variable environment such as natural morphologies perfected by evolution. 

Optimal microstructures of two-phase and multiphase composites are 
drastically different. In contrast with the steady and intuitively expected 
topology of two-material optimal mixture (a strong material always sur- 
rounds weak inclusions), optimal multimaterial structures show the large 
variety of patterns and the optimal topology depends on volume fractions. 
Optimal multiphase structure may contain an enveloping layer but it also has 
"hubs" of a material with intermediate stiffness connected by "pathways" 
(laminate of the best and worst materials) and other configurations that re- 
veal a geometrical esse nce of optimality, see (jCherkaevl . l2009uCherkaev and Zhangl . 



201 ll : ICherkaevl . |2012| ) and Figure [T] Geometries of multimaterial optimal 



structures are not unique, pieces of the same material may occur in different 
places of an optimal structure and they may correspond to different fields 
inside them. Clearly, the method for finding optimal multiphase geometries 
differ from those for optimal two-material structures. 
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Figure 1: Composites made from two isotropic materials and void (white and grey colors 
represent stronger and weaker phases repectively, black corresponds to void) . Comparison 
of geometries of optimal microstructures with maximal effective bulk modulus: (a) two- 
material composite; (b)-(d) three-material microstructures with small, intermediate and 
large volume of stronger material. The amount of weaker material is kept constant in all 
cases. 



In this paper, we follow (ICherkaev and Zhangl 1201 ll ) in developing the 
technique for finding the bounds and optimal structures and we apply it 
to determining optimal elastic composites. The re sults constitute a next 
step from the popular "topology optimization" , see ( lBends0e and Sigmundl . 
20031 ). that is a problem of optimal layout of a material and void. Namely, 



we describe an optimal distribution of two materials and void, or an optimal 
two-material porous composite made from an expensive and strong material, 
a weak and cheap one, and void. Such problem is useful for many multi- 
physics design applications when additional properties besides elasticity are 
important. For example, the presence of one of the materials can be inde- 
pendently required by the reasons different from an optimal elastic response. 
The above-mentioned consideration puts forward new challenges, such as ap- 
proximation of complex optimal structures with simpler ones, selection of the 
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cases where three and more materials should be mixed in the same spot, etc. 



Background of the research Formally, the problem of optimal struc- 
tures can be formulated as a question of minimizers of a variational problem 
with nonquasiconvex multiwell Lagrangians; the wells represent components' 
energies plus their costs and the minimizers (Young measures) are stress fields 
in the materials of an optimal composite. The challenging open problem is 
to build the quasiconvex envelope for Lagrangian with three or more wells. 
The problem is addressed by (i) finding exact bound (the lower bound for 
the quasiconvex envelope) and (ii) approximating these bounds by special 
class of minimizers. By building the lower bound, we also obtain sufficient 
conditions on optimal fields in materials that hin t on the search for geomet 



ric patterns determin ing optimal structures, see (IMiltonl . l2002l : lAlbin et al. 



20071 : ICherkaevl . |2009[ ). 



The existing techniques for th e bound such as Hashin-S htrikman method, 
see (IHashin and Shtrikmaij. 119631): translation method, see ( iLurie and Cherkaevl . 



19821 : iLurie and Cherkaev ^ 



1986 



Tartarl.ll985l:lMiltonl.l2002[) or a nalvtic method 



of Bergman- Milton, see (iBergmanl . Il978l : iMiltonl . Il981al . |2002[ ). produced a 
numbe r of re s ults for two-material mixtures in the last 25 years, see e.g. the 



books flLurid . 119931 : ICherkaevl . l2000l : lAUaird . 120021 : IMiltonl . 120021 : iDacorognal . 



20081 ) for examples. These techniques, however, do not provide exact solu- 
tions for multiwell problems. This, in turn, poses a mathematical challenge 
that needs to be addressed. 

In the last three decade s , the rn ultimaterial optimal composites have been 
studi ed by Milton ( Milton . 1981bh. Lurie fc C herkaev ( Lurie and Cherkaev . 



1985 ). Milton & Kohn ( iMilton and Kohnl . ll988[ ). Gibiansk y fc S i gmun d ( IGibiansky and Sigmundl . 



2000|) among others. In 1995, Nesi published a paper (iNesil . Il995l ) about 



bounds for multimaterial mixtures that are better than Hashin-Shtrikman 
ones. 



In the last years (2009-2012), a new techniqu e was suggested (ICherkaevl . 



20091 : ICherkaev and Zhangl . l2011t ICherkaevl . |2012|) for finding optimal bounds 
for multimaterial mixtures and tested at a couple of examples. The essence of 
the new technique is coupling the translation meth od with the Alessandrini- 
Nesi inequality, see ( lAlessandrini and Nesil . l200ll ). that order and restrain 
values of the fields in any optimal composite. Roughly speaking, in the 
case tackled in this paper (elastic 2D microstructures of maximal stiffness 
for a mixture of two materials and void), it states that the sign of a stress 
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field is constant in the whole microstructure. The technique w as used to 
find t he bo unds of isotropic 2D multimaterial composites, see (jCherkaevl . 



20091 . I2OI2I ): a nisotropic conducting comp osites made from two materials 



and void, see (jCherkaev and Zhang . l201ll ). In all these well as 

in the present research, the lower bound is a multifaceted surface. Its an- 
alytic expression varies depending on the convexity of the wells. In these 
problems, the structure of optimal fields in each material of the composite 
is analyzed and described for all cases but one, where the bound and the 
conjunc ted optimal microstructures are still slightly different, see the discus- 
sion in ( ICherkaev and Zhang] . I2OIII ). Our guess is that there is another, yet 
unaccounted, inequality which becomes active in the latter case. 



2. Problem setting 

2.1. Notation 

Reference to periodic homogenization. Consider a domain Q C 
filled with two linearly elastic materials and a void. Non-homogeneous 
distribution of phases in the domain is determined by its division into three 
disjoint subsets Qi, i = 1,2, 3. Suppose that a boundary value problem 
(BVP) of linearized elasticity is posed in Q. If non-homogeneity of material 
layout is given by a fine partition of the domain then a solution to the BVP 
may be difficult to obtain even numerically. In such case, it is convenient to 
make use of the homogenization theory of periodic media in determining the 
simplified, effective Hooke's law in Q prior to solving the BVP. 

The homogenization t heory defines effective elastic m oduli of a given pe- 
riodic structure, see e.g. ( Cherkaev . 2000l : Tartar . 2000) for the most recent 



exposition of the topic. In this paper, we solve an inverse homogenization 
problem: we find a structure of a multicomponent composite that stores 
minimal stress energy in a given homogeneous stress field. We assume that 
the properties of constituent materials and their volume fractions rrii G [0, 1], 
i = 1,2, 3, mi + 7712 + 1713 = 1, are given. 



Definitions. Due to local character of homogenization, in the sequel 
we consider arbitrary x & Q which is sufficiently distant from the boundary 
dfl. Let Y = [0, 1]^ denote a corresponding unit cell periodically extended 
to R^. Assume that Y is divided into three disjont subcells Yt, i = 1,2,3, 
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whose areas are fixed. Write 



y= [j Yi, m^rrii, J]m, = l (1) 

j=l,2,3 i=l 

and set (61,62) for a Cartesian basis in Y. Let stand for a space of 
plane, second-order symmetric tensors, and for the space of plane Hooke's 
tensors. Next, choose 

1 1 
El = — ^(ei ® ei + 62 ® 62), E2 = —j={e-\ ® 61 - 62 ® 62), 

Y v2 (2) 

-B3 = ^(ei <8) 62 + 62 (8) 61) 
v2 

for the basis in E^. 

Suppose that Y\ and Y2 are filled with elastic isotropic materials whose 
constitutive properties are given by Ki = = l/iii, i = 1,2, where ki 

and /li stand for bulk and shear moduli of i-th phase. Let = — +00 
which means that the third phase corresponds to void. Introduce a set A = 
{Ai, A2, Ayoid} where 

Ai^^Ei^Ei + ^{E20E2 + E3® Es) (3) 

represents Hooke's compliance tensor of i-th non-degenerate isotropic phase. 
In the sequel we assume that the materials are well-ordered, i.e. Ki < K2 < 
K3 = +00 and Li < L2 < — +00. 
Set 

To = 1 61 (g) ei + ^ 62 62 

for the average stress tensor in Y. Components 1 and g denote principal 
values of Tq and (61, 62) stands for its principal basis. It follows that 

To — So El + Dq E2, So — —j^, Do — —jS 

and So I Dq represent spherical and deviatoric components of tq. Without 
loss of generality we may assume < 1. It is understood that tq denotes a 
value (calculated at x e of a stress field solving the homogenized BVP. 
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Next, define a set of stress fields statically admissible in Y 

S= |r: rGL|(y,E2), divr = OinF, J^T{y)dy = To 

where L^(y, E^) stands for the space of L^-functions with values in and 
F-periodic in Q. Fields r = T{y) belonging to this space are thus uniquely 
represented in by one spherical and two deviatoric components, respec- 
tively given by s = s{y) and di = di{y), d2 = d2{y), such that 

Tll+r22 , ru-T22 , rr 

s = — — , di = — — , da = V2 

hence r(y) = s{y) Ei + di{y) E2 + d2{y) E3. 

Due to F-periodicity, r G S is endowed with two properties: 

- function detr(?/), y eY, is quasiaffine hence 

/ det T{y)dy = detro = g; (4) 
Jy 

- function detr(?/), y eY, is locally univalent with detro, that is 

det T{y) > a.e. in F if det tq = ^ > (5) 
and the latter remains valid if " > " is replaced by " < " , 



see (lAlessandrini and Nesi l200ll ). The above-mentioned properties do not 
result in any restrictions on r G S, they simply unveil certain characteristics 
of the stress fields related to assumed F-periodicity. Nevertheless, (jlj) and 
([5]) are of great significance in bounding the stress energy which is the central 
part of the study. 

Decomposing the determinant function of a stress field according to 

2 det r = - {d^ + d^) 

and considering q E [—1,1], allows for rewriting ([5]) in the form 



s\y) > df{y) + d^iy) a.e. in F if G [0, 1], 
s^y) < dfiy) + d^iy) a.e. in F if G [-1, 0]. 



(6) 



For further considerations, let us rephrase the requirements imposed on r G 
S. First, define a set 

^uni = |t: r G L^(F, E^) with univalence property as in ([HD, |. 

Next, write the restriction on the mean stress (/y t = tq) in a form 

Eav = Si, Dij, i,j = 1,2: miSi + 771282 = So, 

rriiDu + m2-Di2 = Dq, miDi2 + m2-D22 = 0, 
S,'>D^, + Dl iige[0, 1], 

S,'<D^, + Dl, iff?G[-l,0] } 

where 

Si = —! s{y)dy, Dij = — [ dj{y)dy, i,j = 1,2, (7) 

denote average spherical and deviatoric stresses in non-degenerate phases. It 
follows that S C Srci where 

Srci = |t: r G Suni and such that Si,Dij G Sav, i,j = 1,2| (8) 

stands for the set of relaxed stress fields, i.e. fields with neglected differential 
constraint divr = in F. 

2.2. Composite materials of minimal stress energy 

Energy bound and extremal effective material properties. The 

(quadrupled) stress energy density in Yi, i = 1,2, is calculated according to 

U.ir) = 4 [r : (A, r)] =K,s^ + U {d? + d^) (9) 

and we set U3{t) = due to assumed r = in void. The contraction r : {Ai r) 
is realized by a standard operation [r]^ (Ai) [t] in the basis ((21). Here [r] and 
[Ai) stand for a vector and matrix representations of respective quantities and 
[r]-^ denotes a transpose of [r]. Effective energy is thus calculated according 
to 



^ i=i 



r G E ^ (10) 



and Uo{g) is bounded from below by 



U.ig) = mi\UoiQ) V-^asinO 



Bounding the stress energy allows for restricting the values of effective 
constitutive properties. Indeed, by introducing K^:, and A^, linked simi- 
larly to ([2]) one may claim ?7*(^) in the form 



U4g) =4 [ro : ( A tq 



1 



-{K4l + gy + L4l-gy 



(11) 

With this notation, K^, and L^, represent coupled bounds on effective moduli 
of a composite for fixed g. They may be understood as constitutive properties 
of a homogenized medium adjusted to the external stress tq = SoEi + Dq E2 
in a sense of storing the minimal amount of energy in two directions Ei, E2 
simultaneously. 

Note that in light of the above, the requirement of isotropy imposed on 
the effective medium is redundant. Indeed, the component of A* related 
to the direction E3 ® E^ may be arbitrary as tq : E^ = 0. Non-isotropic 
microstructres may thus be optimal, i.e. such that the amount of stress energy 
stored in them equals U^{g). Details on this topic are presented in Sec. IH 

Let us find formulae for K^, and L*. To this end, note that by varying 
g G [—1, 1] on the r.h.s. of f lTTj) we obtain a family of functions that are 
quadratic in g and U^{g) represents an envelope of this family. Solving the 
system 

U,{g) - I (k, (1 + + (1 - gf 



d 

dg 



0, 



(1 + gf + (1 - gf 



0, 



allows for determining the coefficients of U^{g). They read 



KM 

L^{g) ■■ 



U*{g) ^ I - g dU^.{g) 



1 + g 



+ 



2(1 + g) dg 
1 + g dU^{g) 



(12) 



1 — g 2(1 — g) dg 



Functions in ( fT2|) are extremal if their values belong to dGmA, i.e. the 
boundary of G- closure of set A. Recall that GmA contains all effective 
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Hooke's tensors obtained by homogenization of components belonging to A, 
taken with arbit rary microstructure and fixed volume fractions rrii, see e.g. 
( ICherkaevl . |2000[ ) for further reference. 



Calculating energy bound by the translation method. In what 
follows we briefly describe a procedure of determining f/*(f)). For this we 
make use of the translation method which proved to be an efficient tool in 



' Cherkaev 




2000; 


Milton 




2002) 



from introducing a translation parameter a e T C M and rephrasing in 
the form 

f/i(r) = Fj(r, a) — 2 a det r, i 



1,2, 



where 



Fi(r, a) = {Ki + a) + (Lj - a) [dl + d^) . 
With dl]) taken into consideration we calculate 

/ Ui{T)dy = y2 F,{T,a)dy-2ga. 

Next, we neglect the differential constraints on the stress field in Y which 
allows for taking the infimum in ( ITOj) on the set enlarged from S to Erei- 
Another effect of dropping div r = in F is that the optimal stress field r e 
Srei can be determined independently in each phase. Moreover, the search 
can be reduced to non-degenarate phases only as r = in void. Consequently, 
one obtains 



^{g,a) = M\^ [ Fi{T,a)dy 



(13) 



r G Srci 

By ([8]) it is possible to split the latter task into two steps. First, we set 



^i{Si, Ai, Di2, a) = inf | J Fi{T, a) dy 



and we continue with 

$(^, a) = min $i + $2 



Si, Dij G Sav ( ■ 



1 = 1,2, (14) 



(15) 
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Stress energy estimation predicted by the translation method is thus given 

by 

U^{g) >Utr{g) = maxf^<^{g,a) -2ga a G T j (16) 

and the equahty U^{g) = f/tr(f?) holds if r G Srei solving ( TT^ . or equivalently 
f|T^ and f|T5l) . proves to be statically admissible, i.e. r G S. This require- 
ment is fulfilled if the stress field components in neighboring materials are 
consistent with divr = in Y. If this is the case then formula for f/tr(^?) is 
optimal, as it corresponds to the boundary of GmA and it may be substituted 
in f[T^ for calculating extremal coupled effective properties of a three-phase 
composite. Explicit calculations of Utr{g) is a subject of Sec. [3l and proving 
its optimality is postponed for Sec. HI 



3. Lower bound on the stress energy: Sufficient condition of opti- 
mality 

We proceed by explicit calculation of r G Srei in two steps defined by (ITll) 
and f|T5|) . This in turn allows for determining Utr{g) by proper adjustment of 
the translation parameter a in f|T6|) . Consequently, bounds on effective con- 
stitutive properties K^{g) and L^{g) related to dGmA are obtained through 
f lT^ under the assumption that energy bound predicted by the translation 
method is optimal, i.e. that f/*(^)) = Utr{g) holds. Discussion of the latter 
is provided in Sections H] and |5l Below we establish the sufficient condi- 
tion of optimality in terms of stress fields related to Utr{g)- With Ki, Li, 
i = 1,2 given, the sought condition turns out to be dependent on mutual 
relations among mi, m2 and g. It results in the division of a polyhedron 
n = {{g, mi, 7712) : g G [—1,1], mi G [0,1 — m2],m2 G [0,1]} into several 
regions of optimality. Table [T] provides a brief guide to the sequence and 
results of calculations and Fig. [2] shows an exemplary cross-section of 11 by 
a plane m2 = const. 
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Table 1: A guide to calculations of the lower bound on the stress energy, related effective 
constitutive properties and sufficient condition of optimality. 



range 




regions 


bound 


KM 


sufficient 




of g 




in Fig. E] 




and L^:{g) 


condition 






1,0] 


A' 


m 


m 


(EH), dMD and 


dBD 






B' 


dMD 




dlZD, dHD and 


dsZD 






C 




dMD 


dlZD, dHD and 


m 






D' 






(dZD, dSl and 


m 




1] 


D 




dSSD'^ 


([HD, d2QD and 


m 






C 






dlZD, dH and 


m 






B 




(m 


dlZD, dH and 


m 






A 


m 


m 


dUD, d^HD and 


m 



unrestricted E 



" For remarks on regions D, D', see Sec. |5l 




-1 -mj 1 Q 



Figure 2: Regions of optimality related to a cross-section of a polyhedron 11 by a plane 
TO2 = 0.35. Constitutive properties of nondegeneratc materials are fixed to Ki ~ 1, 
Li =2,K2^ 3, L2 - 4. 

3.1. Stress fields in phases 

In f fT4D we wish to obtain ^i{Si, Dn, Di2,a) > —00 for i = 1,2, as such 
property is crucial in subsequent derivation of a nontrivial energy bound 
Utr{g). To this end, we ffist discuss the conditions under which the integrand 
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Fj(r, a) is bounded from below by a convex function. Next, we determine 
opti mal relaxed stres s fields r G Sum by making use of the Jensen inequality, 
see (ICherkaevl . I2OOOI . Sec. 1.2). Applied to our case, it states that if Fi(r, a) 
is convex in r then its integral over Yi takes a minimum value on a constant 
stress field being the average of r over yj. Thus, with r decomposed into 
spherical and deviatoric parts, we expect the minimizers to be expressed in 
terms of averages Si, Dn, Di2. 

Assuming that Si, Dij, i,j = 1,2, are prescribed and Ki ^ Li, we consider 
the following cases: 

(I) If a G {—Ki, Li) then all terms in -Fj(r, a) are convex. 
For g E [—1, 1], from the Jensen inequality it follows that 

^i{Si, Dii, Di2, a) = mi{Ki + a) Sf + mi{Li - a) [D^ + D^) . 

Components of optimal relaxed field constant in Yi and 

read 

s{y) = Si, dj{y) = Dij a.e. in Yi. (17) 

(Ila) If a > Li then the s-term in -Fj(r, a) is convex and d-terms are concave. 
For g G [0, 1] we have > + (i| a.e. in Yi hence -Fj(r, a) is bounded 
from below by a convex function. Indeed, 

{Ki + a)s^ + (Li - a) {d^ + d^) > (Ki + Li) 

and the " > " relation above is replaced by " = " when = d^ + (i| 
a.e. in Y^. Consequently, 

MSi, Ai, A2, «) = m,{K,, + L,) (18) 

Components of optimal r G Suni are given by 

s{y) = Si (constant field), ^^^^ 
di{y) = Si cos 6 (y), d2{y) = Si sin 6 (y), a.e. in 

with 6{y) arbitrary in Yi up to the restrictions given in ([7]). 
(lib) \i a = Li then the s-term in -Fj(r, a) is convex and d-terms vanish. 

For g G [0,1] we obtain ^i{Si,Dii,Di2,a) in a form identical to (fTSjl . 
Formulae determining the components of optimal relaxed field read 

s{y) = Si (constant field), ^^^^ 
dj{y) arbitrary up to > df + d^ a.e. in Yi. 
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(Ilia) If a < —Ki then the s-term in Fj(r, a) is concave and (i-terms are 
convex. 

For Q G [—1, 0] we have < + (i| a.e. in Yi hence 

{Ki + a)s^ + {Li - a) [dl + d^) > {K^ + U) {dl + d^) (21) 

and the " > " relation above is replaced by " = " when = dl + (i| 
a.e. in Y^. Consequently, 

^,{Si, Da, A2, a) = m,{Ki + U) [Of^ + A2) • (22) 

Components of optimal relaxed stress field are constant in Yi and given 
by 

s{y) = ±^A'i + Ai dj{y) = D^j a.e. in Y,. (23) 

(Illb) If a = ~Ki then the s-term in -Fj(r, a) vanishes and rf-terms are convex. 
For g G [—1,0] we obtain ^i{Si, Dn, Di2,a) in the form identical to 
fl22]) with optimal relaxed fields 

s{y) arbitrary up to < c/^ + a.e. in Yi, ^^^^ 
dj{y) = Dij (constant field). 

Other relations between g and a are not discussed here as they are irrelevant 
in further study. 

3.2. Calculation of the energy bound: Regions of optimality 

Having $j = ^i{Si, Dn, Di2,a), i = 1,2, explicitly calculated, we now 
turn to the problem of determining ^{g,a) and Utr{g) through (IT^ and 
(ITB]) . Moreover, substituting optimal Si, Dij in the formulae for s{y) and 
dj{y) leads to sufficient optimality condition of the translation energy bound 
in terms of a stress field in i-th phase. 

3.2.1. Case of g E [0,1] and a E {Li, L2) : Regions B and C 

According to the discussion in Sec. 13. Ij items (I) and (Ila), set 

$1 =mi(iri + Li)S2, 

$2 = m2{K2 + a) Si + m2(L2 - a) {D^^ + D^^) . 
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It follows that 

^{g,a)= min{$i + <l>2} 

subject to: mis'! + 171282 = Sq, 

rriiDu + m2-D2i = Dq, 
miDi2 + m2-D22 = 0, 



(25) 



52 > + ^2 



12' 
2 

21 "T -L^22- 



From the KKT optimality conditions we conclude that: (i) D12 = -D22 = 
0, (ii) $1 + $2 is minimized with respect to D21 if -Dii is maximized. All 
constraints in fl25l) are satisfied if we set 

Du = min | —-Do, Si 
[mi 

Consequently, further discussion splits into two subcases corresponding to 
regions of optimality B and C in Fig. [2J Forthcoming results are rather 
straightforward to obtain in both regions, hence we omit the details of cal- 
culations. 

Region B : Assume that 

Dii = —Do < Si. 
mi 

By this, D21 = and (l25l) transforms to 

$(^,a)= mm^^mi{Ki + Li)S^ + m2{K2 + a)Si^ ^^^^ 

subject to: miSi + m2S2 = Sq, 

hence one may replace the above with a one-dimensional unconstrained opti- 
mization problem where 5*1 is treated as a variable. The necessary condition 
of optimality 5($i + ^2)/dSi = allows for determining the function Si{a) 
which is substituted back in fl26l) . and formula for 5*2(0) follows from the 
constraint. The lower bound on stress energy in region B given by 

UtriQ) = max |$(f), a) — 2 ^ a aG(Li,L2)| = 

ii±izM^r?r A' ^^^^ 

= ^— {Ki + Li) + 2gK2 
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results from solving d[^{g, a) — 2 ga]/da = 0. In this way, two critical values 
of a are obtained. The one corresponding to maximum in f l27p reads 

a = - ^ ^ [Ki + Li) - K2. 

2 grrii 

Consequently, optimal mean spherical and deviatoric stress components in 
region B are given by 

_ l + g-2v/g^ _ l-Q _ 

— -E- ; J^ll — —E- 7 ^12 — U, 

52 = v^^^^, ^21 = 0, 1^22 = 0. 

m2 

Sufficient optimality condition of (I27|) expressed in terms of stress fields 
in phases 1 and 2 respectively follow from substituting Si, Dn, D12 in ( IT^ 
and 5*2,-021,-022 in (117]). 

Region B is represented by a curvilinear rectangle P2-P3-P5-P6 in Fig. |2j Its 
boundaries are determined according to the following scheme 

a < L2 =^ mi > 11)^^3(1^2, q), 
a>Li =^ mi < ■il)B^^{m2, q) , 
miSi > Dq =^ Q e [m2, 1] 



where 



^/gm^il + Q- 2^gm2) Ki + Li 
y/grn^{l + g - 2y/grn^) Ki + Li 



(29) 



ipB-Dim2, 0) 



2g K2 + L1 



Region C : Conversely to previous paragraph assume 

Dii = Si< —Do, 
mi 

which results in 

^{g,a)= min |mi(-ft'i + -Li) 5*^ + 

+ m2 [{K2 + a) Si + {L2 - a)Di^] 

subject to: miSi + m2S2 = Sq, 
miSi + m2-D2i = -Dq. 
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The algorithm of calculations is similar to the one presented for region B. It 
follows that the stress energy in region C is bounded from below by 

Utrio) = max a) — 2 f)a; a;G(Li,L2)| = 

(30) 

2 7712 

{Ki + Li)(l - + {K2 + L2)mim2 



with 



2 mi 



a 



{L2 -K2) + ^ [(1 - m2)(Ki + Li) - mi(K2 + L2)] 



2 1^ f)mi 

and 



5*1 — ^= , Dii — 5*1, D12 — 0, 

v2mi 

m2 + g _ m2-g _ 

V2m2 y lm2 



(31) 



By substituting fl3T]) in f lT9|) and f|T7|) respectively we obtain sufficient 
optimality condition of ( l30l) in terms of stress fields in phases 1 and 2. 

Region C is represented in Fig. |2] by a curvilinear triangle P1P2P3. Its 
boundaries are determined by the following expressions 

a < L2 nil > i)A^c{m2, g), 

a>Li rrii < ilJc-E{m2, g), 

TTiiSi < Dq =^ g^[0,m2] 



where 



1712(1 -m2){Ki + Li 

VJA-dm2,g) = 



^c-e("^2, g) 



1712 + g){.K2 + L2) 

m2{l - m2){Ki + Li) 



(32) 



(m2 + g)iK2 + L2) - 2{L2 - Li)g 

3.2.2. Case of g E [0, 1] and a > L2: Region A 

According to the discussion in Sec. 13. H item (Ila), set 

<l>i = mi{Ki + Li) S,^ z = l,2. 
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Hence 

^{g,a) = mill {$1 + $2} 

subject to: rriiSi + 771282 = Sq 

and 

fAr(f?) = max a) — 2 a a>L2^. 

It is immediate that the function to be maximized monotonically decreases 
in a hence we set a = L2. Results obtained in the remainder of this section 
correspond to the region of optimahty A in Fig. [2j 

Region A : Proceeding analogously to previous cases we derive the 
lower estimate of stress energy in region A. It takes the form 

2 mi{K2 + L2) + m2{Ki + Li) 
Optimal values of mean spherical components of stresses read 
1 + ^ K2 + L2 



Si 
S2 



y/2 mi{K2 + L2) + m2{Ki + LiY 

1 + g Ki + Li 

y/2 mi{K2 + L2) + m2{Ki + Li)' 



(34) 



Substituting Si and 5*2 in f|T9|) and fl20l) leads to sufficient optimality 
condition of fl55]) in terms of stress fields in phases 1 and 2 respectively. 
For fixed Ki, Li, i = 1, 2, and arbitrary 777-2, region A is described by 

A = |(^, 777i): < 7711 < i'A-dm2, g) if ge [0,7772], 
< 7771 < iJ/K-B{m2, g) if g e [7772, 1]| 

see fl29|) and fl32|) . Curvilinear sides of A are represented in Fig.[2]by lines P1P3 
and P^Pe- Note that A splits into Ai and A2 with the interface represented 
by a curve P3P7. This division is explained in Sec. 14.2.31 

3.2.3. Case of g e [-1,0] and a G {-K2, -Ki): Regions B' and C 
According to the discussion in Sec. 13. H items (I) and (Ilia), set 

$l= 777i(iri + Li) {Df, + D^,), 

$2 = 7772 (ir2 + a) Si + 7772(^2 - tt) {D^ + D^^) ■ 
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It follows that 



^{g,a)= min{$i + $2} 

subject to: niiSi + m-iS-i = Sq, 

■rriiDu + m2-D2i = Dq 
miDi2 + m2-D22 = 0, 

'21 "T -L^22- 



(35) 



5f < DA + z^A 



From the KKT optimality conditions we conclude that: (i) D12 = D22 = 
0, (ii) 5*2 minimizes $1 + $2 if Si takes its maximal value. All constraints in 
(135|1 are satisfied if we set 

5*1 = min < — 5*0, -Dii 

Similarly to Sec. 13.2.11 the case splits into two subcases. They correspond to 
regions of optimality B' and C in Fig. O 

Region B' : Assume that 

^1 = —^0 < ^11. 
nil 

By this, 5*2 = and f p5|) transforms to 

<l>(^, a) = min |mi(Ki + Li) + m2(i^2 + «) ^2^} 
subject to: niiDu + m2-D2i = -Do- 
Proceeding analogously to previous sections we obtain 

Utr{g) = max a) — 2 f)a; a G (— ^^2, — -^i)| = 

{l-g-2,/^^y (36) 



2 mi 

with 



-{Ki + Li) - 2qL2 



1 v/-gm2(l - g - 2V-gm2) ^^^ 

a = (Ai + Li) + L2. 

2 grrii 
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Optimal mean values of spherical and deviatoric components of stresses 
in region B' are thus given by 



^1 - —r^ , ^11 - ^= , ^12 - U, 

m2 

Sufficient condition of optimality of (1361) result from substituting 5*1, Dn, 
Di2 in ([23D and (HID. 

Region B' is represented by a curvilinear rectangle P'^P'^^Pl^P'^ in Fig. [21 
Its boundaries are determined according to the following scheme 

a > -K2 ^ mi> ipA'-B'{m2, g), 
a < -Ki =^ rrii < ipB'~D'{m2, g), 
miDu>So [-1,-7712] 



where 



y/-gm2{l - g~ 2^/-gm2) Ki + Li 

i^A'-B'{m2, g) = ir~T~r- 

2g K2 + L2 

^/-gm2{l - g - 2y/-gm2) Ki + Li 



(38) 



V'B'-D'("^2, g) 



2g K1 + L2 



Region C : Conversely to previous paragraph assume 

Si = Dii < — 5*0, 
nil 

which results in 

^{g,a)= min|?Tii(Ki + Ll)Dfl + 
+m2[(K2+a)5|+(L2-«)/^|l] 
subject to: rriiDu + 171282 = Sq, 

miDu + 7712-021 = Dq. 

It follows that the stress energy in region C is bounded from below by 
Utr{g) = max |$(^, a) — 2 a a E {—K2, —Ki)^ = 

{K2 + L2)g^ , , 

= —r, '—^{K2-L2)g^ (39) 

2 7712 

^ [Kx + Li)(l - 7712)^ + {K2 + L2)mi7772 

27Hi 
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with 

ry = 

2 



and 



a = ll{L2-K2)-^ [(1 - m2){Ki + L,) - m,{K2 + L2)] ] 



5^ = ^, /^21 = ^, D22 = ^. 

v2m2 v2m2 

Substituting 5*1, Dn, Z^i2 in fl25]) and 5*2, -D215 -D22 in f lTT|) results in 
sufficient optimality condition of (!39|) . 

Region C is represented in Fig. |2] by a curvilinear triangle PiPI^P!^. Its 
boundaries are determined by the following expressions 

a > -K2 ^ mi> tpA'-c{m2, g), 
a < -Ki =^ nil < 'ipc'~E{m2, g), 
miDii < So =^ gE[-m2,0] 



where 



1712(1 -7712) {Ki + L 

i'A'-c'{fn2, 0) = 



ipc'-Eim2,e) 



1712 - q){K2 + L2) 

m2(l - m2){Ki + Li] 



(41) 



(m2 + q){K2 + L2) - 2{L2 + Ki)q 

3.2.4- Case of g E [—1,0] and a < —K2: Region A' 

According to the discussion in Sec. 13. H item (Ilia), set 

$i = mi{Ki + Li) {Dl + Dl) , 1 = 1,2. 

Hence 

^{g,a)= min{$i + $2} 

subject to: niiDu + m2-D2i = Dq, 

miDi2 + ^2-022 = 0, 

and it is immediate that D12 = -D22 = 0. 

The estimate of the stress energy is determined as 



f/tr(^) = max |$(p, a) - 2 a < -K2Y (42) 



Function to be maximized in (|42|) monotonically decreases in a hence we set 
a = —K2. Results obtained in the remainder of this section correspond to 
the region of optimality A' in Fig. [21 
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Region A' : The lower estimate of stress energy in region A' takes the 
form 

Tl ( ^ (1-g)^ {K, + L^){K, + L,) 

UtriQ) = ^ I r ^ I fi^ ^ r \ + 2^^2- 43 



Optimal values of mean spherical fields read 

^ ^ 1-g K2 + L2 

V2 mi{K2 + L2) + m2{Ki + Li) 
l-Q Ki + Li 



(44) 



^/2 mi{K2 + L2) + m2{Ki + Li] 



Considering (jSj) in (|23l) and leads to sufficient optimality condition 
of (H3l) in terms of stress fields in phases 1 and 2 respectively. 
For fixed Ki, Lj, i = 1, 2, region A' is described by 

A' = |(p,mi,m2): < mi < V'a'-C'("^2, g) H g E [-m2,0], 

< mi < ^A'-B'(m2, ^) if ^ e [-l,m2]| 

see ( l38l) and ( HTi) . Curvilinear sides of A' are represented in Fig. [2] by lines 
P1P3 and -Ps-Pg. Note that A' splits into A'^ and Aj with the interface repre- 
sented by a curve -P3-P7. This division is explained in Sec. 14.2.31 

3.2.5. Case of g E [—1, 1] and a G {—Ki,Li): Regions D, D' and E 

The case of a G {—Ki, Li) is discussed in Sec. 13.11 item (I). Both $1 and 
$2 are described by 

= m, [{Ki + a) + (L^ - a) (A'i + AI)] , t = l,2. 

For G [0, 1] it follows that 

«)= min{$i + $2} 

subject to: miS*! + m25'2 = 5*0, 

miDii + m2-D2i = -Do, 
miDi2 + m2-D22 = 0, 
5f > Dl, + Dl„ 
Si > Dl + Dl. 



(45) 
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In case of ^ G [— 1, 0], last two constraints change into 



52 < D^, + Df,, 
Si < Dl + Dl. 



The KKT requirements are that D12 
arbitrary problem (HSj) takes the form 



D 
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0. Hence, in case of 



$(^, a) = min jmi [{Ki + a) S^ + (Li - a) DfJ + 
+ m2 [{K2 + a) Si + (L2 - a) 
subject to: niiSi + m2S2 = Sq, 

rriiDii + m2-D2i = -Do- 
Applying necessary conditions of optimality leads to 

{K2 + a)So ^ {L2-a)Do 



Si 
S2 



mi{K2 + a) + m2{Ki + a) 

{Ki + a) So 
mi{K2 + a) + m2{Ki + a) 



D21 



mi{L2 - a) + m2{Li - a) 

{Li - a) Dp 

mi(L2 - a) + m2{Li - a) 



Technically, determining t/tr(^?) and optimal a G {—Ki.Li) requires sim- 
ilar algorithm to the one used in previous sections. Corresponding formulae 
are not presented here due to their complexity. The limiting cases of a = Li 
or a = —Ki are dealt with separately in the sequel of this Section. As a 
result, three additional regions of optimality appear. Namely, region E cor- 
responds to a G {—Ki^Lij and regions D and D' are related to a = Li and 
a = —Ki respectively. 

Region E : The problem is now to find the formulae for two curves 
which separate region E from D and D'. These curves are given by implicit 
functions '?/'d-e('^i, ''^2, = and V'd'-eI^i, ''^25 = where 



V'd-e("^i,"^2, q) ■ 
^D'-Eimi,m2, g) 



d_ 
da 
d_ 

da - 



a) — 2 ga 
$(^, a) — 2 ga 



a=Li 



a=-Ki 
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2 , 2^(1 - mi) 



It follows that 

V'd-e("^i, g) = ^D-E{mi, m2)(l + + 

nil 

V'D-E(mi,m2) = 

_ m2(Ki + Li) [(mi + m2)(ii:i + Li) + 2 mi (7^2 - i^i 

2mi [mi(is:2 + ^i) + "i2(i^i + ^i)]^ 

and 

ipD'~E{mi, m2, ^) = ^d'-e("^i, "^2)(1 - e')^ + ^-^^^ 

mi 

^D'-e("^1,"^2) = 

m2{Ki + Li) [(mi + m2){Ki + Li) + 2mi(L2 - 



^ ^ 2 



2 mi [mi (fsTi + L2) + m2 (iCi + L 
Region E is bounded by the hnes 
mi = 1 — m2, 

mi = ipc'-E li ge [-1,0], 

mi = ipc-E if ^1 e [0, 1], 

^d'-e("^i,"^2, g) = 0, 

V'd-e("^i,"^2, g) = 0, 

see (l32l) and (HT]) . Curvilinear sides of E described above are respectively 
represented in Fig. [2] by lines P1P3, P1P3, Pl^P'^ and P2 Pi- 
Region D : Assuming a = Li and repeating the discussion for region 
E or region B (Sec. 13.2. ip we conclude that the stress energy estimation and 
optimal mean fields in D are given by 

rr + {Ki + Li){K2 + Li) 

^"^'^ - mi(ir2 + Li)+m2(iri + Li) - ^^^^ 

and 

{K2 + Li){l + g) l-g _ 

i~>i — — T, ^11 — — ;= , ^12 — u, 

V2[mi(A'2 + Li) + m,(A'i + Li)] V2m, 

V2 [mi(is:2 + ^1) + ^2(7^1 + L{ 
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Region D is represented in Fig. [2] by an area bounded by two straight 
lines: (i) g = I, (ii) mi = 1 — m2 and two curves: (i) '?/'d-e("^15 '"^25 = 
(line P2P4), (ii) rm = ^b-d (line P2P,). 



Region D' : Assuming a = —Ki and repeating the discussion for region 
E or region B' (Sec. I3.2.3P we conclude that the stress energy estimation and 
optimal mean fields in D' are given by 

(1-g)^ iK, + L,)iK, + L,) 
2 mi{Ki + L2) + m2{Ki + Li) 

and 

S^ = l±l. c^^- {K, + L2){1-q) ^ 

^ y2mi' " ^/2[mi{Ki + L2) + m2{Ki + 



S2 ^ 0, D2. ^ (^1 + ^0(1-^) ^, 



(49) 



Region D' is represented in Fig. [2] by an area bounded by two straight 
lines: (i) q = —1, (ii) mi = 1 — m2 and two curves: (i) ^/'d'-e('^i, ^^2, f?) = 
(line P^'^i), (ii) mi = t^b'-d' (hue P^P^). 

Sufficient condition of optimality of fj46|) result from substituting fH9l) in 
fl20|) and ([17]). Similarly, considering (gg]) in ([21]) and ([I?]) leads to sufficient 
optimality condition of ( l48l) . 

5". 5. Bounds on effective isotropic properties 

Making use of (fT2l) allows for calculating bounds on effective isotropic 
properties in each optimality region where U^,{g) is determined. From the re- 
sults obtained in the preceeding section and by assuming that f/tr(f?) = f^*(f?) 
it follows that formulae for K^{g) and L^{g) can be derived in any region ex- 
cept E. It also has to be pointed out that the equality of energy formulae is 
only conjectured in regions D and D'. Consequently, -R'=^(^?) and L^:{q) corre- 
spond to dGmA (boundary of the G-closure of the set of effective properties) 
only in the high-porosity regions A, B, C and A', B', C. Recall that K^,{q) 
and -L*(f?) are related to dGmA only if optimal stress fields predicted in Sec. 
13.11 are statically admissible. In Sec. ID we prove that this is the case for 
high-porosity regions. We conjecture the same property for D and D', see the 
discussion in Sec. [51 
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Region A : 



mi ^ 1712 



Ki + Li K2 + L2 



- L2, L^{g) = L2 



(50) 



Region A' : 

K^g) = K2, L,(q) = 



nil _|_ 



Ki + Li K2 + U 



K2. (51) 



Region B : 



K,{g)^K2- 



[(1 + g)^/grn2 - 2 ^] [l + ^ - 2^ q 7112 



2 mi ^(1 + q) 



{Ki + Li), 



T(\ v/g^(l + g - 2^/gm^) 

= ^ (-^1 + Li) - K2. 

2 nil Q 

Region B' : 

\/-Qm2{l - Q- 2^-Qm2) 



(52) 



KM 



L*(^) = L2 + 



[(1 - g)y/-gm2 + 2 ^] [l - ^ - 2y/-gm2 



(53) 



2 mi ^(1 — g) 



(Ki + Li). 



Note that ^ < in B'. 
Regions C and C : 



L*{g) = 7; 



{K2 - L2 



{l-m2Y 
mi(l + g) 



{Ki + Li) 



1712 + 



{K2 + L2) 



(1 — m2)^ ^ , r \ , ^2 ~ 



mi(l - g) 



{Ki + Li) + 



m2(l - £») 



m2(l + g) 

{K2 + L2) - {K2 - L2) 



(54) 



Region D : 

K.{g) : 



mi 



m2 



Ki + Li K2 + Li 



Li, L,{g) = Li. (55) 
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Region D' : 

KM ^ Ku LM ^ + ^) " - A-, (66) 



i A:^(e), L«(Q) 

^HS 




\ 

\ 

\ 




Figure 3: Comparison of optimal bounds K^,{g) (solid lines), Lf{g) (dashed lines) and the 
Hashin-Shtrikman estimates ifnsi -^hs (dotted lines). Values of functions are calculated 
for mi = 0.17, m2 = 0.35 and Ki = 1, Li = 2, K2 — 3, L2 = 4. Symbols qac, Qbd, Qa'C 
Qb'D' refer to the anisotropy level of tq at the interfaces between respective regions; g = 
at the interface between A and A', g ~ 777,2 at the interface between C and B, g ^ —7772 at 
the interface between C and B'. 

Figure [3] illustrates the comparison of functions K^{g) and L^{g) repre- 
senting coupled lower bounds on isotropic properties of a three-phase com- 
posite in different regions with the Hashin-Shtrikman uncoupled bounds 

Li, 

2Ki + Li. 



Kns 



nil 



+ 



1712 



Ki + aK K2 + ax 



-1 



ax, ax 



mi 



+ 



m2 



Li + ai L2 + ai 



-1 
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Estimates K^s, Lhs are independent of g E [—1, 1], as they do not incorporate 
an information on the anisotropy of Tq. 

Note that K^:{g) < K^ls for all g E [—1, 1] and K^:{g) = K^is in region D 
while the inequahty < Lhs is slack in all regions, see the discussion in 

Sec. EH 



4. Optimal high-rank laminates in high-porosity regions 

In this Section we show that optimal relaxed stress fields determined in 
regions A, B, C and A', B', C coincide with statically admissible stress fields 
r G S. The task is two-fold. First, we make use of the differential constraint 
divr = in deriving additional requirements on r G Efgi. Next, we show that 
these requirements are fulfilled in certain microstructures, so-called laminates 
of high rank. 

4.1. Compatibility of stresses on phase interfaces and mean stress fields in 
rank-one laminates 

In calculations of optimal r G Srd, the differential constraint divr = 
in Y (equilibrium equation) is neglected. Consequently, energy-minimizing 
stress fields are determined in each phase independently. It follows that 
components of optimal relaxed fields may be incompatible with divr = on 
material interfaces which in turn means that r ^ S. 

Suppose that two materials meet in a given microstructure at a line F 
and let n and t denote a normal and tangent to F. In the sequel we consider 
microstructures where phases are arranged in layers hence F takes a form 
of a straight line. Moreover, we assume that stress field in each layer is 
constant. By this we claim that if a given non-degenerate phase Yi, i = 1,2, is 
distributed in p layers Yi^i, Fj 2, • • ^i,p, then optimal r is layer-wise constant 
in Yi. It follows that if p = 1 then r is constant in entire Y^. Equlibrium 
equation is thus fulfilled identically in each phase. 

Constraint divr = requires that [[Tn]]r = where [[-jjr denotes a jump 
on the interface F between layers of materials. If we set r^, m = 1,2, for 
constant fields on both sides of F then the jump condition may be rewritten 
in a form 

(n - T2) : (n ® n) = {n - T2) : (n^t) = (ri - T2) : (t ® n) = 0. 
It follows that stress fields with (ri — r2) : (t^t) 7^ are con ipatible with the 



equlibrium constraint hence statically admissible in Y, see ( ICherkaevl . 12000 
Sec. 14.2.2) for full discussion of this topic. 
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Figure 4: Rank-one connectivity of stress fields Tm, m = 1, 2, 3 such that : (ei (8)62) = 0: 
(a) simple laminate L(12) with phases taken in proportions f3i, 1 — /3i, a normal to the 
interface F given by n = ei and stress fields ti and T2; (b) simple laminate -^(13) with 
phases taken in proportions P2, ^ — P2, a normal to the interface F given by n = 62 and 
constant stress fields ti and ra; (c) graphical interpretation of compatibility conditions 
(ti — T2) : (ei ei) — 0, (ti — T3) : (e2 (81 62) — and mean fields ^^(12) in laminate i(12), 
T^(i3) in laminate i(f3). Vectors Ei, E2 are defined in ([2]). 

Here we discuss stresses ri, T2 in two materials arranged in a rank-one 
laminate -^^(12). Compatibility of stress fields in 1/(12) is also referred to 
as rank-one connectivity at T. Let ri and T2 denote rank-one connected 
stress fields in materials layered in proportions (3 and I — (3 respectively. 
Resulting mean field takes a value Ti(i2) = (3ti + {1 — /3)t2- Examples of 
rank-one connected stress fields and their mean values in simple laminates 
are sketched in Fig. |H High-rank laminates are constructed by repeated 
rank-one layering scheme under the assumption that the materials resulting 
from previous laminations are homogeneous. These type of structures are 
considered in the subsequent section. 



29 



4-2. Stress fields in micro structures 
4-2.1. Regions C and C 

Region C : Continuing the discussion in Sec. 13.2. H one may notice that 
the assumption 5*1 = Du enforces 9{y) = a.e. in Yi. Consequently, optimal 
stress field r G Srei is constant in Yi. We thus calculate 



r2 



m2 



mi 
1 ei (g) ei + 



Q 

m2 



a.e. in Yi, 
62 eg) 62 a.e. in Y2. 



(57) 



In order to prove statical admissibility of ( !57l) we check the compatibility of 
stresses on material interfaces in a L(13, 2) laminate, see Fig. |5l 




Figure 5: Optimal microstructure in regions C and C: (a) layout of materials; (b) layering 
scheme leading to tq = 7^(13^2) for region C (a mirror image with respect to ei (8) ei results 
in the scheme for region C). Stress fields in layers of strong and weak materials are 
represented by white and grey squares respectively, r = in void is represented by black 
square. Circles represent stress fields in laminates. 



Outline of the layering scheme: 
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1. Substructure L{13) is formed: phase 1 and void are laminated with 
rii = 62 and volume fractions Pi, 1 — (3i respectively. Homogenized 
stress field in L(13) is given by t"l(i3) = /3i ri. 

2. Final structure L{13, 2) is formed: phase 2 and L{13) are laminated 
with n2 = ei and volume fractions (32, ^ — P2 respectively. Fields Ti{i3) 
and T2 are rank-one connected if (t/^(i3) — T2) : = 0. Stress field 
in the final structure tl(i3^2) = P2 ^2 + (1-/^2) t"l(i3) satisfies ri(i3,2) = tq- 

Parameters of optimal laminate: Compatibility conditions reduce to 



.,„„:(e,».0 = l ^ A = — 6 10,1], ^^^j 

tl(13,2) : (e2 ® 62) = ^1 =^ (32 = m2 e [0, 1] 

and it is immediate that the constraints on volume fractions of phases in Y 
given by 

(3i{l - (32) = mi, /32 = 1712 

are satisfied identically. 



Region C : In addition to the considerations in Sec. 13.2.31 and due to 
assumed Si = Du we set s{y) = di{y) = Si a.e. in Yi. Further discus- 
sion reduces to the one presented above with g E [—1, 0] taken into account. 
Spherical and deviatoric components of mean stress in each phase are given 
by the same formulae in both regions C and C, see fl?Il) and fHUl) . Con- 
sequently, stress field in laminate L(13,2) fulfills the sufficient condition of 
optimality also in region C with phase volume fractions given by ( 158|) . Layout 
of materials and scheme of layering corresponding to region C are sketched 
in Fig. O 

4.2.2. Regions B and B' 

Region B : Here we continue the discussion in Sec. 13.2.11 with the as- 
sumption Du < Si. Function 9{y), y G Yi, may vary in Yi hence the stress 
field in material 1 are rank-one connected with zero stress in void if d{y) = 
or 9{y) = 71. Taking this into consideration we subdivide phase 1 into two 
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layers, i.e. we set Yi = Yi i + Yi 2- Formulae for stresses read 

1 + g — 2 Wgm2 

Til = Ci Ci a.e. in Yi i, 

mi 

1 + g-2^gm2 . , . 

Ti 2 = 62 ® 62 a.e. m Yio, (,59; 

nil 



J gni2 

T2 = (6i ® 6i + 62 ® 62) a.e. in Y2. 

m2 

Next, we prove statical admissibility of (jSH]) by checking rank-one con- 
nectivity of stress fields in a L(13i,2,132) laminate, see Fig. O 




Figure 6: Optimal microstructure in regions B and B': (a) layout of materials; (b) layering 
scheme leading to tq = ''■l(i3i,2,139) for region B (a mirror image with respect to ei (8> ei 
results in the scheme for region B'). Stress fields in layers of strong and weak materials 
are represented by white and grey squares respectively, r = in void is represented by 
black square. Circles represent stress fields in laminates. 



Outline of the layering scheme: 

1. Substructures L(13i) and L(132) are formed. In L(13i), the first layer 
of phase 1 (field ri^i) and void are laminated with ni = 62 and volume 
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fractions 1 — /3i respectively. Stress field in L(13i) is given by 
^1,(131) = In ^(182), the second layer of phase 1 (field ri 2) and 

void are laminated with n2 = ei and volume fractions P2, ^ — ^2- Stress 
field in L(132) reads tl{i^^) = [^2 ti,2- 

2. Substructure L(13i,2) is formed: phase 2 and -L(13i) are laminated 
with 77-3 = ei and volume fractions /Sa, 1 — (3^. Fields Ti{i3i) T2 are 
rank-one connected if (tl(i3j) — T2) : (ci ci) = 0. Stress field in the 
substructure is given by rL(i3j,2) = /^s ^^2 + (1 - /^s) t"l{13i)- 

3. Final structure Iv(13i, 2, 132) is formed: laminates L(132) and L(13i, 2) 
are layered with ^4 = 62 and volume fractions /34 and 1 — /54. Rank-one 
connection between stress fields holds if (ri(i32)—ri(i3^ 2)) : {.^2®^2) = 0. 
Formulae rL(i3j,2,i32) = P^^oa + (1 - /34)^o,3 and TLi^vz^, 2,1^2) = ^0 link 
the fields in substructures with the effective stress tensor. 



get the following 



(^L(i3i) - T2) : (ei (g) ei) 


= 


» /3i 


(^L(132) - Tl(13i,2)) : (e2 


® 62) = = 


» /32 


^L(13i,2,132) : (e2 ® 62) = 


^' 


» P3 


•7"L(i3i, 2,132) : (ei ® ei) = 


1 





Parameters of optimal laminate: Collecting the combatibility conditions we 

2^1 ' 

(33 S2 

2 Si ' 
V2g 

S2 ' 

^2 

Parameters /3j G [0, 1], i = 1, ... ,4, as a/2 < S'2 < 5'i in entire region B. 
Indeed, it is a matter of straightforward calculations to check that the first 
inequality immediately follows due to ^ > m2 in B and the second one reduces 

to 

mi < il)B-D[Q^m2) 

Ki + Li 

by substituting relevant formulae from ( l28l) and ( l29l) . 

Constraints on volume fractions of phases in Y are satisfied if + 
1^1,2! = mi and 1^2! = m2 or, equivalently, 

(l-/33)(l-/34)/3i + /34/32 = mi , 
(1 - /34)/33 = m2 
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from which we have 



Pi — — rr~ 7; — N 5 P2 



m2(l + g -1^Qvn2) ' 1 + ^ - 2 ' ^^^^ 

/Jg = ^gm2 , (3i = l- J— . 



Note that //(13i, 2, 182) (optimal in region B) morphs into -L(13, 2) (optimal 
in C) at the boundary between regions. This is concluded by setting g = 1712 
in . . . , above. 

Region B' : From the discussion in Sec. 13.2.31 and the assumption that 
< Dii it follows that phase 1 is subdivided into two layers such that 

Yi = Fi 1 + Fi 2 with s{y) = Du a.e. in Yi i and s{y) = —Du a.e. in Fi 2. In 

this way, stress fields are given by 

1 - g-2y/-gm2 . 
Til = Ci (g) ei a.e. m Yi 1, 

mi 

1 - g -2^y-gm2 . 

T12 = 62 ® 62 a.e. m Ki 2, 

mi 

J-gm2 , , . 

T2 = (ei ® ei - 62 ® 62) a.e. m r2- 

m2 

Thus described stress field in laminate L(13i, 2,132) is statically admis- 
sible in region B'. Phase volume fractions are given by fl60|) with g replaced 
with —g. The details of calculations are omitted here as they follow the 
same pattern as presented above. Layout of materials and scheme of layering 
corresponding to region B' are sketched in Fig. [6l 

4.2.3. Regions A and A' 

Region A splits into two subregions Ai, A2 with different optimal mi- 
crostructures. 

Subregion Ai : We use the results of Sec. l3.2.2] to prove that the stress 
fields 

Ti = V2 Si ci (g) ci a.e. in Yi, 

T2,i = a/2 5'2 ei (g) ei a.e. in ^2,1, 

5*2 ~l~ fl S2 — fl 

T2,2 = — ^ ei ® ei H ® ^2 
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are statically admissible in laminate L(123, 2), see Fig. [71 Phase 2 is thus 
subdivided into layers F2,i and ^2,2- In the sequel we assume r2,2 : (ei(8>ei) = 1 
from which it follows that /i = — 82- Values of Si and S2 are given by 




Figure 7: Optimal microstructure in regions Ai and A'^ : (a) layout of materials; (b) layering 
scheme leading to tq = Ti(i32,2) for region Ai (a mirror image with respect to ei iSi ei 
results in the scheme for region A'j^). Stress fields in layers of strong and weak materials 
are represented by white and grey squares respectively, r = in void is represented by 
black square. Circles represent stress fields in laminates. 

Outline of the layering scheme: 

1. Substructure L(123) is formed in two steps: (i) phase 1 and void are 
layered with rii = 62 and volume fractions /3i and 1 — /3i (in this way 
1/(13) with ro,i = /3i ri is obtained), (ii) first layer of phase 2 and i^(13) 
are laminated in the same direction with volume fractions f32 and 1— /32; 
this leads to L(123) with rj^(i23) = 1^2 ^2,1 + (1 - /32) Pi n- 

2. Final structure L(123, 2) is formed: second layer of phase 2 (field r2,2) 
and L(123) are laminated with = ei and volume fractions and 
1 — /Sa- Rank-one connectivity requirement reads (ri{i23) — ^2,2) '■ (ei ® 
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ei) = 0. Stress field ri(i23,2) = /Ss T2,2 + (1 - /33)tl(i23) matches the 
effective tensor by 7-^(123,2) = ''"o- 

Parameters of optimal laminate: Rank-one connectivity requirements sim- 
plify to 

(rL(i23) - 1-2,2) : (ei ® ei) = = (T^^^ ~ 5^ 

7-L(i23,2) : (e2 (^62) = g ^ P3 ^ 



V2S2-I 

and constraints on volume fractions in Y are fulfilled if l^il = mi and |1^2,i| + 
|>^2,2| = m2, i.e. 

(1 - /32)(1 - /33)/3i = mi, (1 - P,)P2 + /33 = ^2. 

Therefore, optimal laminate is parameterized by 



^ _ mi ^ _ m2(l - m2)K - mi(m2 + ^i) 
l-^^2' {1 + g)[{l - 1712) K - nil] 

^ _ g{mi +m2K,) ^ _ Ki + Li 
{1 + g - m2)K, - mi' K2 + L2 



(61) 



Note that k e [0, 1]. 

Our next claim is that /3j G [0,1], i = 1,2,3. Indeed, f3i coincides with 
that in region C and the remaining conditions can be reduced to 

(32 >0 if mi < V^A-c("^2, g), 

/32 < 1 and /53 > if mi < (1 + ^ - m2)K, (62) 
< 1 if mi < (1 — m2)/t. 

For the definition of ipA-c{m2, g) and ipA-B{m2, g) (used below) see Sec. 13.2.11 
It follows that the first constraint in (|62l) is most restrictive. Hence, A = 
Ai U A2 where 

Ai = |(mi, g) : < mi < V^a-c("^2, g), g e [0, 1]|, 

A2 = |(mi, g) : ^A-c("^2, g) <mi< V'a-b("^2, g), g e [m2, 1]|. 

Laminate L(123,2) (optimal in Ai) morphs into L(13,2) (optimal in C) 
at the boundary between regions. This follows from substituting mi = 

V^A-c("^2, g) in (32 and (^3. 
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Subregion A2 : We make use of "the coating principle" , see ( lAlbin et al. 



20071 . Th. 9), in determining optimal microstructure in subregion A2. Lam- 
inate L(13i,2, 182) (optimal in region B) is coated with a layer of phase 2, 
in the direction = ei normal to the interface, and volume fractions 1 — P5 
and /Ss- In this way, L(13i, 2, 182, 2) is obtained. Phases 1 and 2 are thus 
subdivided according to Yi = Yi^i + Yi_2 and Y2 = l2,i + ^2,2 respectively. 
Formulae for stress fields in phases read 

Ti^i = \/2Si ei ® ei a.e. in Fi^i, 

Ti,2 = a/2 5'i 62 62 a.e. in Fi,2, 

r2,i = ^(61 (g) 61 + 62 (g) 62) a.e. in 1^2,1, ^^^^ 

V2 

5*2 ~l~ f2 S') — ^2 

r2,2 = — ^ ei ® 61 + ^-^=—62 ® 62 a.e. m ^2,2- 

Technically, sufficient condition of optimality in subregion A2 do not re- 
strict stress field in 1^2,1 to be spherical. However, the assumed form of 
T2,i proves to be optimal as it is shown in the sequel. We also assume 
■72,2 : (61 ® 61) = 1 which gives /2 = ^2 - S'2 in 
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Figure 8: Optimal microstructure in regions A2 and Aj: (a) layout of materials; (b) layering 
scheme leading to tq = Ti(i3i,2,i32,2) for region A2 (a mirror image with respect to ei (g) ei 
results in the scheme for region Aj). Stress fields in layers of strong and weak materials 
are represented by white and grey squares respectively, t = in void is represented by 
black square. Circles represent stress fields in laminates. 

Outline of the layering scheme: 

1. Substructure L(13i,2,132)is formed along the scheme discussed in case 
of region B hence it is not repeated here. 

2. Final structure L(13i, 2, 132, 2) is formed: second layer of phase 2 (field 
^2,2) and L(13i, 2, 132) arc laminated with = Ci and volume fractions 
1— /Qs and (3^. Rank-one connectivity condition is given by (Ti(i3i,2.i32)~ 
T2,2) : (ei ® ei) = 0. Stress field in L(13i, 2, 132, 2) is linked to the 
effective tensor by rL(i3i ,2,132,2) = p5r2,2 + {i- p5)ro,4 and ri(i3^,2,i32,2) = 

To- 

Parameters of optimal laminate: Rank-one compatibility of stress fields in 
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substructures and resulting relations between laminate parameters read 

(tl(i3i) - 7-2,1) : (ei O ei) = = ^ , 

P3S2 

(tl(132) - ^L(13i,2)) : (62 ® 63) = => /^s = (64) 

(tl(i3i ,2,132) - "^2,2) : (ei (8) ei) = =^ ^4 = 1- — 

and the constraints on volume fractions are given by |yi,i| + |yi,2| = ttii and 
I ^2,1 1 + 1^2,2! = m2, or explicitly 

(1 - ^5) [(1 - Ml - ^3) A + = mi, 

(1-^5)(1- W3 + /55 = m2. 

Introducing 

X — K {1 + g) > 0, y = 2 (mi + m2 k) > 0, x > y (65) 

gives 



_ X [(m2 + ^1)1/ - 2 m2 x] 
^3 = 



(x - yy (1 - /^s) 



m2x'^ - gy^ 

J5 - 



(66) 



(x - y)2 

and we proceed to show that e [0, 1], i = 1, . . . , 5. 

Prom K e [0, 1] it follows that /3i fulfills the condition. Parameter falls 
into required interval if S2 > a/2 in subrcgion A2. This condition is satisfied 
if 2 mi < {1 + g — 2 1712) n. Observe that 

, . {l + g-2m2)K 
mi < ipA-B{m2, g) < 

holds for m2 < g. Thus (3^ e [0, 1]. 

We can assert that /^s > if 1712 x'^ — gy"^ > 0. This inequality can be 
reformulated to mi < ipA-B{m2, g) and the assertion follows. In order to 
prove (3^ < 1 we show that 

{x -yf -m2x'^ + gy"^ >0 (67) 
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in entire A2. To this end, we first rewrite f l67p in the form 



[l + ()){y - yi ){y - y2) > 0, 
yi^2 = (1 ± a/(1 + g)m2 - g) k 



yi > y2, 



(68) 



and the proof falls naturally into two parts. If (1 + g)m2 — g < then the 
roots in (!68l) do not exist and (!67l) follows immediately. 

Conversely, let us assume that (1 + g)m2 — g > 0. Next, make use of (!65il 
and (1681) to calculate 



where my represents values of mi corresponding to yi. To show that ( 167|1 
holds it is sufficient to check if my(m2, g) < 'ipA-c{fn2, g)- The latter is fulfilled 
if (1 + g)m2 — g > —Tn2 which may be concluded from the assumption. 

Comparing the expressions in (j66ll we deduce that the discussion of /^s > 
may be reduced to proving (m2 + g)y — 2m2 x > 0. It is straightforward to 
compute that this inequality is equivalent to mi > ^/'a_c(^2 5 g)- For checking 
if /^s < 1 we write it in a form (x — y) [x (1 + 7712) — y {1 + g)~\ > 0. From 
fl65|) we see that it suffices to show that x {1 + m2) — y {I + g) > 0. This 
requirement reduces to 2 mi < (1 — m2)K which is valid in entire A2 due to 



The property of (32 G [0, 1] follows from f]^ G [0, 1], see (^^. 

Recall that a boundary between regions Ai and A2 is given by mi = 
ipA-c{m2, g), g > m2. By substituting this formula in (3i, . . . , (3^ one may 
check that -L(13i, 2, 182, 2) (optimal in A2) morphs into L(123, 2) (optimal in 
Ai). By the same token, setting mi = '?/'a-b("^25 in /3i, . . . , /Ss leads to the 
conclusion that L(13i, 2, 182, 2) smoothly changes into L(13i, 2, 182) (optimal 
in B). 

Similarly to A, region A' also splits into Ai' and A2'. Optimal microstruc- 
tures in both subregions are the same as in Ai and A2 respectively. 



my(m2, g) 



1 - 2 m2 + + g)m2- g 
2 



K 



mi < ^/'a-b("^2, g) < 



(1 - m2) K 



2 
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Subregion Ai' : For proving optimality of the laminate L(123, 2) we 
make use of the results obtained in Sec. 13.2.41 Stress fields in phases read 

Ti = Dii ci ® ei a.e. in Yi, 

T2,i = -D21 ei ® ei a.e. in 1^2,1, 

9i + D21 gi - D21 . 
T2,2 = 1= — ei (g) ei H -= — 62 ® 62 a.e. m ^2,2 

V 2 V 2 

where Du, D21 are given by (jSj). We assume that r2,2 : (61 ® 61) = 1 from 
which it follows that gi = \/2 — D21 ■ Further calculations are similar to those 
presented for subregion Ai. They lead to formulae for optimal lamination 
parameters written in ( 16T|1 with q replaced by —q. 

Subregion A2' : Laminate -L(13i, 2, 132, 2) proves to be optimal also in 
subregion A2'. Stress fields are given by 

Ti^i = ^2 -Dii 61 ® 61 a.e. in Yi^i, 

n,2 = v^-Dii 62 ® 62 a.e. in Fi_2, 

r2,i = — F^(ei (g) 61 - 62 ® 62) a.e. in ¥2^1, 
v2 

g2 + D21 g2 ~ D21 

T2,2 = 1= — ei (g) 61 H -= — 62 ® 62 a.e. m Y2,2 

V 2 V 2 

and we assume that r2,2 : (ei ® 61) = 1 hence g2 = \/2 — D2i- Optimal 
lamination parameters are derived similarly to those in region A2. They read 

A = ^, & = ^, 

^ X' [{1712 - Q)y' -21712 x'] n _ 

^ m2 (xQ' + 

where x' = (1 — f?) k and y' = y. 

Proof of /3j G [0, 1], 2 = 1, . . . , 5 in subregion A2' follows the pattern set 
in the discussion regarding A2. 



5. Final remarks 

5.1. Remarks on low-porosity regions D, D' and E 

The detailed description of the mentioned regions D, D' and E of large 
volume fraction of the first material, or, equivalently, low-porosity regions, 
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will be provided in a separate paper. Here we restrict ourselves with some 
brief remarks outlining the current results. 



Region D: The Hashin-Shtrikman bound on bulk modulus At 

the boundary of regions B and C that correspond to maximal allowed vol- 
ume fraction mi, the optimal translation parameter reaches the value of Li. 
The energy bound f/tr(f?) in region B transforms into the classical transla- 
tion bound which corresponds to the Hashin - Shtrikman bound o n the bulk 
modulus fo r isotropic composites. This bound is realizable, see (ICherkaevl . 



20091 . |2012[ ). The anisotropic translation bound is attained on certain mi- 



crostructures only when the anisot r opy l evel is not too large, compare the 
discussion in (jCherkaev and Zhangl . 120111 ). The optimal structures for both 
conducting and elastic composites are similar, they are determined by high- 
rank orthogonal laminates L(13i, 2, 132, 1, !)• These structures are obtained 
by enveloping the nucleus laminate L(13i,2, 132) - optimal for the region B 
- by two orthogonal layers of the first material. It is shown in (lAlbin et al. 



20071 ) that such enveloping is stable with respect to the translation bound: 



if the nucleus satisfies this bound, then the enveloped nucleus also satisfies 
it. 



Region D': Similarly, at the boundary of regions B' and C, the optimal 
translation parameter reaches the value of —Ki. However, in this case, the 
energy bound f/tr(f?) in region B' does not give rise to the Hashin - Shtrikman 
bound on the shear modulus for isotropic composites. 

Indeed, Utr measures the energy of a composite subjected to an arbitrary 
stress field whose anisotropy is controlled by ^ G [—1, 1]. Consequently, if 
we set ^ = 1 then the effective energy is optimized only in a direction of the 
applied field tq = [(1 + g)/2] Ei which is spherical, i.e. isotropic. 

On the contrary, setting g = —1 does not lead to a similar conclusion 
because applying the deviatoric field tq = [(1 — p)/2] E2 and retaining the 
isotropy of a composite medium by controlling its response in the direction 
i?3 at the same time is impossible. 



Region E: Guessed optimal structures Optimal L(13i, 2, 132, 1, 1) 
structures degenerate into the most anisotropic L(13,2, 1), when the ani- 
sotropy of the external field increases (the value of g decreases from 1 to- 
wards 0). When the anisotropy level increases even further, the translation 
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bound is not realizable by the known structures. Moreover, it is definitely 
not optimal for strongly ani sotropic structures; the reasons are discussed in 



(jCherkaev and Zhang . 1201 ll ) 



We conjecture that the region E of large volume fractions mi and strongly 
anisotropic loadings correspond to the limiting structures L(13,2,l). The 
bound for this region is presently unknown, and we guess that it corresponds 
to another inequality that becomes an equality in that region. To support 
our guess we mention that: 

- the optimal structure in Region C is L{13, 2). The L{13, 2, 1) structures 
degenerate into them, when the fraction of external layer of the first 
material vanishes; 

- the best known bounds for extremely anisotropic structures [g 

0) correspond to the same str ucture L{13, 2, 1) in that region, see 



(ICherkaev and Gibianskyl . Il996l ): 



- the structures that realize the translation bound for moderately aniso- 
tropic loadings, also degenerate into L(13,2, 1); 

- the L{13, 2, 1) structures degenerate into ^^(13) when the fraction of 
the second material disappears; 

- the L{13, 2, 1) structures degenerate into L{12) when the fraction of 
the third material (void) disappears. 

In the absence of the bound, one cannot prove the sufficient optimality 
conditions for the guessed structures and therefore the global character of 
their optimality. It can be numerically shown, however, that the relative 
gap between a rough bound for the energy and a structure of this class is 



very small, see (jCherkaev and Zhand . l201ll ). Therefore, these structures are 



either optimal or a close approximation of optimal, and can be treated as 
optimal for practical purposes. 
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